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Abstract. The van dcr Waals and Casimir-Polder interaction energy of an 
atom with an infinitely thin sphere with finite conductivity is investigated in 
the framework of the hydrodynamic approach at finite temperature. This 
configuration models the real interaction of an atom with fullerene. The Lifshitz 
approach is used to find the free energy. We find the explicit expression for the 
free energy and perform the analysis of it for i) high and low temperatures, ii) 
large radii of sphere and ii) short separation between an atom and sphere. At 
low temperatures the thermal part of the free energy approaches zero as forth 
£^ , power of the temperature while for high temperatures it is proportional to the 

Cw ■ first degree of the temperature. The entropy of this system is positive for small 

3 ' radii of sphere and it becomes negative at low temperatures and for large radii of 

^^ the sphere. 
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1. Introduction 



Van der Waals dispersion forces plays an important role in different physical, biological 
as well as chemical phenomena Q] [2j O [4] . These forces arise between objects due to 
^ correlated quantum fluctuations. They depend on the shape of the body and their 

structure. The last achievements in Casimir effect and Casimir-Polder interaction have 
been discussed in great depth in books and reviews (see, for example recent book [3]). 
Of special interest is dispersion forces between microparticles and macroscopic bodies. 
In the case of interaction between particle and plate it is commonly referred to as the 
Casimir-Polder force [5]. At short range the energy rises as third power of inverse 
distance between microparticle and the plate. The retardation of the interaction 
should be taken into account at large distances and the interaction energy falls down 
as the fourth power of distance. At separations larger than a few micrometers thermal 
effects become dominating. 

Cavity quantum electrodynamic effects inside a spherical cavity has been 
considered in Refs. [SE1IE]- The non-retarded van der Waals potential for an atom 
inside and outside a metallic bubble was considered in Ref. 9 . The corrections 
for proximity force approximation of the Casimir interactions of an object inside 
and outside a spherical metal shell were calculated in Ref. [TU1 HJ. The perfectly 
conductive sphere with an atom was investigated in Ref. [12]. In Refs. [13l H"4] 
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the Casimir-Polder interaction between an atom and a magnetodielectric sphere was 
considered on the basis of point-scattering techniques. The consideration was related 
to results obtained in Ref. [TS] for interaction an atom with curved surface and Refs. 
[161 117] for an atom with dielectric sphere. The Casimir and Casimir-Polder energy 
were considered in the Refs. p~8l [19] in the framework of the plasma (hydrodynamic) 
model of the sphere which was developed in the Barton's papers [20] [21] in the context 
of the fullerene molecule Ceo- The Barton's approach was based on the Fetter paper 
[22] which was devoted to the properties of the two-dimensional free-electron gas. 
This model was exploited in Ref. [53] for flat sheet. Complete list of references about 
particle in the spherical geometry may be found in reviews |24i 125] . 

The spherical plasma shell was considered in the framework of this model in 
Ref. [IB] in context of the Casimir energy and in Ref. [H] for the Casimir-Polder 
interaction of an atom with spherical plasma surface. In the context of this model 
the conductive sphere is represented by the two-dimensional surface with free-electron 
gas on this surface. The dynamic of this gas is described by hydrodynamics. All 
information about electron gas is encoded solely by parameter with dimension of wave 
number: il = Airne 2 /mc 2 , where n is surface density of electron gas, e and m are the 
respectively charge and mass of electron. It was shown that the interaction energy 
between an atom and the sphere has the form of the Casimir-Polder energy for an 
atom and a plate if atom is situated close to the sphere. Away from the sphere the 
energy has the form of the atom-atom interaction energy with an expression for the 
polarizability of the sphere. The numerical simulations were made for the sphere with 
parameters of the fullerene molecule Ceo and hydrogen atom. 

In the present paper we consider the thermal Casimir-Polder interaction energy 
for a system comprising of an atom and conductive sphere. The sphere is described 
by the plasma (hydrodynamical) model which was considered in detail in the Refs. 
[20l [21] and we use the Lifshitz approach [26] [27] to calculate the free energy of this 
system. We outline the Lifshitz approach in the Sec. [2] 

The calculations performed in the paper could hnd a potential application to 
interaction microparticles with molecule of fullerene. The Casimir-Polder interaction 
of the fullerene molecules C^o and Cm with gold and silicon nitride surfaces was 
considered in detail in Ref. [28] . The fullerene molecule is graphene shaped as a sphere 
and the more realistic model to calculate the Casimir-Polder interaction between an 
atom and a fullerene is the Dirac model. In this model the frequency spectrum 
has linear form and it is valid for low energies. Some calculations related with the 
Casimir and Casimir-Polder energies were made in the framework of this model in the 
Refs. [29] [30j [31] . The frequency spectrum in the hydrodynamic (plasma) model has 
quadratic form and the model is high-frequency approximation of the theory. Both 
models are in qualitative agreement with respect to the Casimir-Polder interaction. 
The calculations performed in the Ref. [31] showed that the Casimir-Polder force 
calculated in framework the Dirac model is smaller then for hydrodynamic model. 

The paper is organized as follows. In Sec. [2] we discuss the Lifshitz approach and 
derive the expression for the free energy and represent it in different forms. Section [3] 
is devoted to consideration of the specific limit cases. We consider the ideal case with 
infinite radius of the sphere but finite distance between the surface and the atom and 
obtain the first correction term over inverse radius of the sphere. We obtain also the 
low and high temperature expansions and consider the case of the short distance from 
the sphere. The expression for the entropy is found in Sec. [4] and it is analyzed in the 
limits of small and high temperatures. We show plots of the numerical calculations for 
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the system comprised the hydrogen atom and sphere with parameters of the fullerene 
Cqq . In Sec. [5] we discuss the results obtained. 

2. Lifshitz approach for an atom near the plasma spherical shell 



We adopt here the approach developed by Lifshitz in Refs. 26, 27 . This approach was 
applied to calculation the Casimir-Polder interaction energy of an atom and plasma 
sphere at zero temperature in Ref. [19]. Let us shortly discuss the approach. We 
consider a conductive infinitely thin sphere with radius R inside the vacuum spherical 
cavity with radius L = R + d which is into the dielectric media with parameters fx, e 
(see Fig. [I}. Due to the spherical symmetry of the problem under consideration the 
electromagnetic field may be divided on the two polarizations which usually called 
as TE and TM modes. We have two concentric spheres and we should consider the 
boundary conditions on the two spherical boundaries. The full information about 
the problem is encoded in these boundary conditions. The information about plasma 



e, |ll^1 ,.<0f~^y^ s, |J=1 




© 



Figure 1. Wc adopt the Lifshitz approach. I: a - the plasma sphere of radius _R, 
b - a spherical cavity of radius L = R+d inside the infinite material with e, fi ^ 1. 
We rarefy the media e = 1 + iizNa + 0(N 2 ) and obtain configuration II: plasma 
sphere with single atom c at the same distance d from the sphere as was cavity's 
surface. 

sphere is encoded in single parameter Q — Aime 2 / mc 2 with dimension of the wave 
number, where e is charge of the electron, m - its mass and n is the surface density 
of electrons on the sphere. The applications of this model for vacuum fluctuations of 
field see in Refs. [20l EU O Q2] . 

The spectrum of the electromagnetic oscillations, u) — kc, in this configuration 
may be found from these boundary conditions (see Ref. [19] for more details): 

Ste = - i Ih'(z £ )$te - -^=H{z e )^ TE \ = 0, (la) 

H(z e )V TM --^H'(z e )* TM \=0, (16) 



J TM 



^2 



where 



^te(^) = J(z) + -J(x)[J(x)Y(z) - J{z)Y{x)}, (lc) 



:r 



*tmW = J(z) + ®J\x)[J'{x)Y(z) - J(z)Y'{x)}, (Id) 
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and J(x) = xji(x), Y(x) = xyi(x), H(x) = xh\ (x) are the Riccati-Bessel functions. 
Here x = kR, Q — flR and z e — Z\fe, z = kL = uiL/c. 

The free energy, J 7 , of the system is expressed in terms of these frequencies in the 
following form 3 : 



1=1 n=l 



In 2sinh — -^4 + In 2 sinh ■ 



2k B T \ 2k B T 



(2) 



where index n numbers the solutions of the boundary condition equations ([2]). Then 
we convert the sum over n to the contour integral 



J- = fc B T^(2Z + l)— j> lnf2sinh^^jd(lnSTM + lnS T E), 



(3) 



where the contour 7 counter-clockwise encloses all real positive solutions of the Eq. 
(|2|). To put the contour to the imaginary axes we should clear up the zero frequency 
behaviour of functions S. If Q 7^ we obtain (y = I + 1/2) 

S TE |^0=U + |;U-^, (4a) 

v , n Kl + !)[(! + 0(1 - £ 2 ) + tt 1+2 '(Z + £ 2 (Z + l))] !_« _i±H .... 

Stm w ->o=Q 7^ X £ 2 , (40) 

where \ = L/R = 1 + d/i?, Q = QR and £ is dielectric permittivity. In the case Q = 
we have 

EteU-s-o = £ 5_ , (5a) 

E T m|u)-x) = ^ -^ X e ^ . (56) 

Therefore, no poles appear at the point w = and we may shift the contour 7 to the 
imaginary axes. 

In the imaginary axes the function under logarithm has zeros at the Matsubara 
frequencies 

.2wnk B T 

w„ = it, n = * z j n = 0, ±1, ±2 . . .. (6) 

a 

We put the contour to the right of the imaginary axes, 

F = k B TY{2l + \)— In 2sinh— — d(lnS T M+lnS T E), (7) 

^ 2tti 7 e+ioo V 2fc B 77 

and then take the limit e — > 0. Integrating by part we obtain 

r=-k B Tj2(2l + l)-l + ^ — coth— (In^ln^)^. (8) 

Then we shift the contour 7 to the imaginary axes. We obtain infinite sum of integrals 
in the imaginary axes between the Matsubara frequencies and integrals over semi- 
spheres around these frequencies. The integrals between £„ and £, n +i and between 
£_„ and £_ n _i cancel each other and we obtain (u = I + 1/2) 

00 +00 
F = k B T^ E ( lnS TM+ln£ T EL =l?M . (9) 

l — l 71 — — OQ 
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To calculate the free energy per unit atom, J 7 ", let us rarefy the media with 
e(iui) = 1 + 4:TrNa(iui) + 0(N 2 ), where a is polarizability of the atom and the density 
of the dielectric matter N — > (see Fig. [T]). In this case the free energy, J", is expressed 
in terms the energy per unit atom T a by relation 

/>oo 

F = Nl F a ^{R + r) 2 dr + 0{N 2 ). (10) 

Jd 



hi 
From this expression it follows that 



d d T 

N~a 4nN(R + df 



^—li-TrW^W (ii) 



After straightforward calculations we obtain 



_fcs7_^ ±25 f kajiuj) ka{iuj) \ 

~( R +Q 2 k \hS G ™ g tm )^ • (12) 



= 1 n— — oo 

where 

x Jte{ik) 



2 -\ Q S ' 2 (x)e> 2 (z) + s' 2 (x)e 2 (z)^P 



G tIi =: - e 'i( z ) s 'i( z ) ~ ei(z)si(z 
and 



/ 7TX I *2*X 

si(x) = d-^Ii+i/2(x), ei(x) = \j—K l+l/2 {x) (13) 

are the modified Riccatti-Bessel spherical functions. 

For renormalization we throw away terms which survive in the limit Q — > and 
arrive with formula 

OG +00 



2k B T£l ^ v ^ / 



(R + d) 2 



s 2 



(x)e 2 (z) s' 2 {x)e' 2 {z) + s' 2 {x)e 2 {z) v - 



fTE{ix) JTM{ix) 



(14) 



where x = tuR/c and z = luL/c = w(R + d)/c, 

/te(«x) =H si{x)ei{x), 

x 

fTM{ix) = 1 sj(x)ej(a;), 

x 

are the Jost functions on the imaginary axes, and the sum with prime means that 

zero term has factor 1/2. The formula (|T4"|) has the same form as for zero temperature 

case obtained in Ref. [19] but with summation over Matsubara frequencies instead of 

integration. 

Now we use single-oscillatory approximation for a: 

a{iw) = -5^-5. (15) 
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At Matsubara frequencies we have 

ft \ gg X a (°) na\ 

a(«U) = ~^j2 = 72 • (16) 

^i % + 1 % + 1 

In the problem under consideration the three different effective temperatures appear 

hui a he he 

" = 2^' R = 2ri^R> d = 2^d' {7} 

The variables are expressed in the terms of these temperatures by relations 

£„ 2irk B Tn T 

n—, (18a) 



0J„ 



fw a T w ' 



£„i? 27rk B TRn T 

x = = = n—, (18o) 



c 



he Tr 



£ n L 2irk B TLn ( T T\ , x 

Let us estimate these temperatures for molecule Cqo and hydrogen atom. For molecule 
C-eo EH] we have R = 0.342nr7i and the frequency w Q for hydrogen atom reads 
[5^ 1 155 | IMj u„ = 11.65eV = 17.698 • 10 15 Hz. Taking into account these parameters 
we obtain T w = 2.15 • 10 4 K, T R = 1.06 • lO 6 ^. 

Now we obtain another representation of the formula (flU . We rewrite Eq. (TTi| 
in the following form 

V ' J=l n=0 



where 



W = -f^, (196) 

c' 2 I ^\ n' 2 ( -,\ I r,' 2 !' 



, 2 (x) e 2 (z) <(x)er(z)+ s r(^K(^)^^ riQ , 

JTE(IX) JTM{IX) 

Here a = T^/T and r = <i/i?. The function g\ consists of two parts which comes from 
TE and TM polarizations of electromagnetic field, correspondingly. Hereafter we will 
denote these contributions as gJ E and g™ . 

The function gi(n) has no poles in the imaginary axes. Taking into account the 
Abel-Plana formula we obtain 

V"Vz?f ^ r°° 9i{t)dt 7T gi(ia) + gi(-ia) 

, ^lyn) = 



f^, Jo t 2 + a 2 2a e 2 ™ -1 

dt. (20) 



n=0 

a + t 



i f°° (9i(it) - gi(-it)\' 



2a J n V e 2rrt - 1 / , 



a — t 



In the first term we change integrand variable t = ak/(cw a ) and the contribution of 
this term is exactly the same as was obtained in Ref. |19] for zero temperature 

he f°° 

Eo = pan Z ^2 / dka(iuj)e (k), (21) 

itR 6 {l + ry J Q 
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where 

oo 

eo(*) = q£) "#(*)• (22) 



/=i 



The second and third terms give the temperature corrections 

hw a 



T _ a(0)ei(a) 

•>! — 



■CJ- 



i? 3 (l + r) 2 e 2 ™ 
a(0)huj a 



1 



7ri? 3 (l + r) 2 J Q 
a(0)huj a 



62 ft) 

3 27Tt ^ 

e 2 (a£) 

-,27TCtt 



In 



a + t 



In 



a — t 

l + t 



where 



ei(t) 



7ri? 3 (l + r) 
-Q ^ i/ [a(*t) + «(-«)] = 5Re (it) 



1-t 



dt 
dt, 



(23) 



;=i 



&2{t) = —Q^2^[gi{it) - gi(-it)} = $Xe (it). 



i=i 

We note that the functions gi(±ia) and gi(±iat) (also ei(a) and e2(ai)) have no 
dependence on the temperature. Indeed, at the beginning (J2J they depend on n via 
a; = nT/TR and z = uT/Tr(1 + r), only. After applying the Abel-Plana formula 
we should replace n — >• ±ia (or n — > ±iai). Therefore, we obtain that x — > iq a 
and z — > ig a (l + r) and whole dependence on the temperature is accumulated in the 
exponential factor the Eq. (j2"31 via the single parameter a — T u /T. 



3. Limit cases 

3.1. An atom near the ideal plate: 51 — > oo and R — > oo 

First of all we have to compare our formulas with well-known case of an atom near 
the flat surface. There is a need to take the limit of infinite conductivity 51 — > oo 
to obtain ideal surface case and then take the limit of infinite radius of the sphere, 
providing finite distance, d, between an atom and the surface. To take limits we use 
the expression for free energy in the form of Eq. (|14j) . Taking the limit of infinite 
conductivity 51 — > oo we obtain the following expression for the free energy 

2k B l sr^ ' 



T a = - 



R? X 2 ^ 



'"KJ^i 



VX 

1=1 

„/2, -, ,2/ 



s 2 (x)e 2 (z) 3 ' 2 {x)e' 2 {z) + s' 2 {x)e 2 {zY- 
si(x)ei(x) 



s'i{x)e[{x) 



(24) 



ui=e„ 



where x = ^ n R/c, z — x\ and ;^ = l+r = l + d/R. 

To obtain expansion over r = d/R <C 1 we use the uniform Debye expansions for 
modified Bessel functions from textbook |35j . Taking into account these expansions 
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X n=0 1=1 y I 

X J 1 + lb (3[f(2;) ~ t{Z)] + [t3(X) ' t3{Z)] + ^^^ + ' 



(25) 



where t(x) = 1/vl + x 2 and ?j(a;) = l/t(x) + In y+jtjtjt- Now we replace summation 
over v = 1+1/2 to integration over v and change the integrand variable v to fc_L = zV-R. 
We observe that the uniform expansion over \jv corresponds to the expansion over 
1/R. For this reason, in the second term of expansion we may set z = x up to first 
power of 1/R and the expression in braces takes the form 1 + 1 5 (£„/ ck±) /2k±R. We 
expand the exponent in above formula over d/R <C 1 and obtain 



e-^HW] = e - 2 ^V 1+ ^T |i + £*_ t ( Jj*_ j + . . . J . (26) 

Taking into account these expansions we obtain the following expansion the free 
energy up to R ; 

i2? r°° ( d h A d 2 k 2 *) 

T a = -2k B T Y, '<*<#*) / dk ± k ± q n e~ 2d ^ 1 - 3 | + -^ + ^ , (27) 

^ Jo { R 2Rq?i RQn J 

where q£ = fc^ + £ 2 /c 2 . The first term of this expansion corresponds to the case of an 
atom near flat ideal surface and it has the same form as in Ref. |36j . The next terms 
give corrections due to finite radius of the sphere. 

Let us discuss an error which appears when we replace summation by integration. 
The precision of calculations is down to R . The second term of expansion over \jv 
in Eq. (|25p gives contribution ~ R^ 1 and when we replace summation by integration 
the error will be outside the considered precision. Let us consider the first term in Eq. 
125): 



Yv^v 2 + z2 ^ 2uWz) ~ v(x)] - (28) 

2 = 1 

To estimate error we use the Abel-Plana formula: 

E /(") - -n\) + f /m* - * jf m l:J}T ) dt (29) 

with f(y) — v\Jv 2 + i2 e -2t/[?7(2:)-?7(x)] > j^. j g eaS y ^ snow that the term /(|) gives 
contribution ~ R~ 2 which is outside of the precision and we drop this term. The 
domain of the integration in the last term is to be divided for three intervals namely, 
0<t<x,x<t<z = xx and t > z. Considering each interval separately we obtain 



2tVz 2 - t 2 e- 2 ^ , f 2Wz 2 - t 2 e~ 2 ^ cos(2v 2 ) , f°° 2Wt 2 - z 2 sin(2v 3 ) , 
-dt+ I „_, : 2> dt+ I e2fft + / 3) dt, 

(30) 
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where the functions read 

fi(t) = v z 2 — t 2 — \J x 2 — t 2 — t\ arctan . — arctan ■ 



/a(t) = \J z 2 — t 2 — t arctan 



V Vx 2 - 1 2 Vz 2 - 1 2 

Vz 2 - t 2 



v 2 (t) = tln ^t 2 -. 



t'3 



(t) = y/t 2 -z 2 - \ / t 2 -x 2 +tln 



zt + Vt 2 



xt + Vt' 



2 _ -,2 



Now we take the limit R — > oo and expand above integrals over 1/ R. The expansion 
of the first integral reads: xe~ T /24, where t = AnkBalT /he and x = £ n R/c. Because 
of the factor i? -3 in the free energy (P5| the contribution of this integral is ~ R~ 2 . 
Both rest integrals in (|30|) give exponentially small contributions ~ e _27rx . Therefore, 
the error which appears when we replace summation to integration is order R~ 2 and 
it is outside the precision under consideration. 

3.2. Low temperature expansion 

Let us consider now the low temperature expansion: T <C T u , Tr, and T <C Tp_ + d- The 
last relation means that we consider distances r = 4 <C t#. The first temperature 
correction given by Eq. (j2"5)l gives exponentially small contribution, T^ ~ e -^ T ^/ T _ 
To find expansion of the second temperature correction (|23[) we take into account that 
the function e2(t) depends on the tT/Tn and tT/Tn(l + r). Owing to the exponential 
factor we have to take into account only small domain of integration t < 1 in the 
integral. Let us find the power expansion of the function giiit) over the t and keep 
the main terms. 

1. TE contribution. The expansion begins with s 5 : 

e TE (t)= - s s Q ^yg±QM + ..., (31) 

6(3 + Q)V 

where s = tT/T R and \ = 1 + r = 1 + d/R. 

2. TM contribution. The expansion begins with s 3 : 

™ frt _ 3 2 5 15Q + 72(5 + Q)x 2 +40Qx 4 + 88Q X 7 
62 W "" S 7 60Q? + -"- (32) 

This expression holds for T/Tr <C Q- The case Q = should be considered separately. 
We set Q = at the beginning and obtain 

„tm, _nf,3 44 .107+184* 



=0-01^46-^-^^)^-0. (33) 

Therefore, the main contribution for the case of the finite conductivity, Q > 0, 
comes from TM mode and it reads 



jT = 2a(o)hu a r 3 /-°° / t 3 ■ h] 



tt X 6 R 3 T^Jo \e M ~l 



a + t 



a — t 



dt. (34) 
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Because of a — T^/T > 1 we expand logarithm over t <C a and obtain 



Therefore, 



and 



where 



t> 



1 



a + t 



a — t 



t* 



2t , 
— dt 



1 



1 / 1 a 



120a 






Hca{0) T 4 _ 4tt 3 a(0) (k B T) i 



In dl ■■ ■■ I I — 1 —dl : : - — — • (35) 

(36) 
(37) 

ST = 45(1 + r)6' (38) 

and Ecp = —3hca(0)/8ird 4 is the Casimir-Polder energy. The plots of Sq and St are 
shown in the Fig. [5] for the case of the molecule Cqq and hydrogen atom. 



607ri? 4 (l + r) 6 T 4 15 (1 + r) 6 {he 



J^ a — Ecp { Sq + St [ — J 



2r 4 



0.005 



0.004 



0.003 



0.002 



0.001 



0.000 




8 S 



Figure 2. The plot of the functions Sq and St from Eq. I|37ll calculated for 
molecule Ceo and hydrogen atom. The function Sq is the zero temperature 
contribution and depends on the parameters of the sphere and the atom. The 
function St is universal and depends on distance r, only. 



3.3. High temperature expansion 

In this case it is better to use Eq. ([14 



J • = - — — 2^^ 2^ u{i£n) 



(R + dr^ 



n=0 



s?(x)ef(*) , ^We^W + ^We^,); 



fTE(ix) 



fTM(ix) 



, (39) 



w=£„ 



where x = fci? and z = fcL = k(R + d), and sum with prime means that zero term has 
factor 1/2. 

At high temperatures the Matsubara frequencies £„ ~ nT —¥ oo for n > and 
these terms in sum over n give exponentially small contribution ~ e~ 2 ^" r . The main 
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contribution comes from zero mode, n = 0, and it maybe calculated in manifest form. 
Indeed, taking the limit n —> in above expression we obtain the following contribution 

^ a _ k B Ta(0) " v{l + 1) 6r 4 + 24r 3 + 33r 2 + 18r + 4 

x2i? 3 Z^ x 2(i+i) KB J ttW 2r3(r + l) 4 (r + 2)3i?3 ' W 

where x = 1 + d/R and r = d/i?. We have the following asymptotic for short and 

large separations between an atom and the sphere 

_ k B Ta(0) d_ 

4d 3 ' R^ ' 

3k B Ta(0)R 3 d 

J- = — — — S> 1 

d 6 ' i?^ 

Extracting the Casimir-Polder energy we obtain 

T a = E CP S T ^, (41) 

J-R 

where 

_ _ 2r(6r 4 + 24r 3 + 33r 2 + 18r + 4) 

T = 3(r + l)4(r + 2)3 ' (42) 

5.^. Short distance behaviour, r = 4 < 1. 

Let us consider Eq. (f2"3"|). In the limit r — > the series is divergent for great I. For 
this reason we use standard formulas (9.3.1) from Ref. |35) : 



_/„(*)«* (££)", (43a) 

Y v {z) «--=(-) , (436) 



>2lTV 

2 

which are fulfil for i^ — > oo. Taking into account these formulas in the expression for 
gi(ia) we obtain 



^' ^ "' l }] (4^ 2 -l)Q-8^2 ) ^ Q = 0,r<a 2 



for TM contribution, and 



2^ + Q ~1 -2^_!l, g^ TO , 



^{^ E M+ft TE Ha)}«-^T7T«<i ^-' ^." (44) 

Q 



for TE contribution. We make summation for finite Q > and r < Q and obtain 

TM 1 - 2 X 2 + 9 X 4 ^ 1 . , 

61 (Q) " 4 X 2 ( X 2 -1)3 "4^ (45a) 

e TE (fl) w ^| (1 _ xarctanh I) « Q 2 lnr (456) 

2 X 4 

For the case of the ideal surface, Q — > oo, the TM contribution will be the same, but 
contribution of TE mode changes from logarithm to first inverse power 
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Therefore, for short distances d <C R and finite conductivity Q > (or Q — > oo) 
the main contribution comes from TM polarization and it reads 



a(0) huj a 

e 1 



*r—^-&-. («) 



If Q = we have 



e7%) « 1 ~ 5x2 ~ 17x4 ~ 27x ° « -- g- (47) 



and therefore for Q — > we obtain 



^r 3Qa(0)i? fia; a 



as should be the case. 

Let us consider Eq. (j2"31 . By using the same asymptotic as above we observe that 
the series over I in function e 2 (a£) is convergent for any r including r = 0. Indeed, in 
the limit I — > oo we have 

- iv {gj^iat) - gj^-iat)} » &f (f^ ( 4 9) 

- iv { g ™ { iat)-gl™(-iat)} * - ^(g^. (50) 

Therefore the series is convergent as Z - ' even for r = (x = 1) and this term gives 
finite contribution and we may threw away it comparing with terms considered above. 
The term with zero temperature contribution (1211) maybe considered in the same 
way. In the limit r — > the series is divergent for large I. For this reason we use the 
Debye uniform expansion formulas from Ref. |35) in which we change z — > zjv and 
make the limit v — > oo: 

w *> "vKr(s)'- (51 °» 



Taking into account these formulas in the expression for gi(k) we obtain 

3-1-2*/ 

"*<*> ^ ^ + 2W (52) 

Now we integrate over k with factor Q in Eq. (|2~Tj) 

dk 1 n \J~Q 



(53) 






, k 2 + k 2 a pQ + 2(kR) 2 2k a ^Q v + V2^q a 
and obtain series 

,3 ,,-l-2f 
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in which we take limit of short distances, x = 1 + r — > 1. If Q ^ we obtain 

y* * v \,-l-*v _ 7T 1 ~ 2X 2 + 9X 4 7T 

tl^ka X 8k a X 2 (X 2 -1) 3 8fc a r3- [ ° 0) 

Therefore, we arrive with expression 

fto; a a(0) 
^° " 8rf^' (56) 

which is valid for r — d/R <C Q = ili?. In the case of Q = we obtain from Eq. (fMf 
that E Q = 0. 

Therefore, for Q > and r = d/fi < Q we obtain free energy 

^ = -^i^" + ^rir (57) 

This expression is in agreement with high and low temperature expansions. 
Indeed, taking the corresponding limits in (|57[) we obtain 

r a =- k -^,T->oc, (58a) 

^=-^,T->0, (586) 



which coincides with Eqs. (|4T|) and (|37|) . 

4. The entropy 

Taking the minus derivative with respect to the temperature (with ks) from the free 
energy given by Eqs. ([23| we obtain expression for entropy S = S\ -\- S% {a = T u /T): 

5 ag) eiCo) ( ™ A 2 (59fl) 

i? 3 (1 + r) 2 VsinhTra/ v ; 

Q a(0) 1 

J2 — 



r°° i a. t / 7ra \ 2 

,, / {te; + e 2 (l-27raicoth7rat)}ln — dt. (596) 

R A TT(l + r) z J 1—t Vsmh7ro/ 

The entropy at low and high temperatures maybe found from Eqs. (|36l) and (|40p 
namely, 

16ir 3 q(0) /fc a T\ 3 2 Q (0)fcg /T\ 3 
5= ~ (TT^b^J = 15 (1+7)* U; J ' T«T fl) T w ,T R+d , (60a) 

q(0) 6r 4 + 24r 3 + 33r 2 + 18r + 4 
5= ~R3- 2r3(r + l)4(r + 2)3 ' T>T R ,T u ,T Il ^. (606) 

The typical behaviour of the entropy is shown in the Fig. [3] for sphere with 
parameters of the fullerene Cqo and hydrogen atom with distance between them 
d = R/2, (r — 0.5). For this situation we have the following numbers: R = 
3.42i, Q = QR = 4.94 ■ 10~ 2 , a(0) = 0.667i 3 , q a = uj a R/c = 0.0202. 

From the Fig. [3] we observe that the Casimir-Polder entropy tends to zero as third 
power of the temperature and it is zero at zero temperature according the Nernst heat 
theorem. It is positive for arbitrary temperatures. In Ref. [35] was observed the 
negative sign of the entropy for short separations, r — > 0, and low temperatures in the 
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SxIO 10 




0.10 T 



Figure 3. The plot of the entropy S calculated for the hydrogen atom and 
the sphere with parameters of the fullerene Ceo arl d distance from fullerene 
d = R/2 (r = 0.5). The entropy for large interval of temperature is shown in Fig. 
a). In Fig. b) we show the entropy for low temperature. The exact numerical 
calculation is upper curve and lower curve is approximate expression (160 £>l ) - 



case of the atom with plate. We note that the low temperature dependence given by 

hc/2-KkBR- To compare with the case of atom near 



Eq. (|60a|) is valid for T <^T R 

the flat boundary considered in Ref. [36] we should take the limit R — > oo and in this 

case the relation T <Tr= hc/2irkBR no longer be valid. 

To show appearance the positive entropy for finite radius of sphere and negative 
sign of entropy for infinite radius of sphere we use expansion over 1/R obtained above 
(|27p . For the case of static polarizability, a(iu) = a(0), we obtain from Eq. (l27l) 



where 

'/o 



l« 



+ 



2re T 



(e T 



T a = e cp [ m + -gm 



r 2 e T (e r + 1) 
(e T - l) 3 



(61) 



Vi 



-: 1 



T 
16 



-1)2 

2re r 
(e- - l) 2 4 
x dt coshi 



iT 2 e 



2 c rl 



2(e T - l) 3 
r 3 (r 2 -4) 



r 3 e r (e 2T +4e T + l) 

2(e T - l) 4 
dt coshi 



i T /2 t sinlTi 48 J t/2 t sinh^i 

and t = 2T/Td = AnkBdT/hc. The function rj was found in Ref. [36], the function 
771 gives the correction due to finite radius of sphere. Now we use these expressions to 
calculate the entropy. Taking the minus derivative with respect to the ksT we obtain 
the following expression for entropy 



3q(0) 

s = ^dF a > 



(62) 



where a = r\' Q T + rr/[ T . The numerical calculations of the entropy for large R by 
using formula (JHTJ) for different values of r = d/R = 0, 0.05, 0.1 are shown in the Fig. 
0] Numerically we found that for r = d/R > 0.085 the entropy becomes completely 
positive. 
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Figure 4. The plot of the function cr calculated for different values of r = 
d/R = 0, 0.05, 0.1. The plot a) shows the function for large domain of 
r = inkgdT/hc. In Fig. b) we show the function cr close to the origin. Starting 
with r = d/R = 0.085 the entropy becomes completely positive. 



5. Conclusion 



In above sections we considered in detail the thermal Casimir-Polder interaction of 
an atom with spherical conductive surface which models fullerene. For sphere we 
use plasma (hydrodynamic) model [201 121] : the sphere is two-dimensional conductive 
surface with finite conductivity. All information about sphere is encoded in one 
parameter Q, = Aitne 2 / mc 2 with dimension of wave number. To obtain the free energy 
we adopt the Lifshitz approach in framework of which the spherical shell is situated 
inside the spherical cavity into the infinite dielectric media. Rarefying this media we 
obtain the free energy per unit atom. Renormalization procedure is simple - we throw 
away all terms which survive in the limit O — >• because this case corresponds to 
the case without sphere and the energy should be zero. The opposite case Q — > oo 
corresponds to the ideal case. The general expression for the free energy is given by 
Eq. dfU) or Eq. ([25]l. 

We made analysis of this expression for different cases. In the limit of infinite 
radius, R, of the sphere with finite distance, d, between an atom and sphere we obtain 
the expected expression for the energy for the case of the atom close to the plate. We 
also found the next term of expansion ~ d/R. In the problem under consideration we 
have three different parameters with dimension of the temperature: 



1L 



27rfcs 



T r , 



he 



2nk B R 



T d 



he 



2nkBd 



(63) 



Each temperature connects with specific scale namely, first connects with scale of atom, 
the second - with scale of sphere and the third connects with distance between an 
atom and the sphere. For molecule of fullerene Cqo we have T^ = 2.15 ■ 10 4 K, Tr = 
1.06 • 10 6 K. For low temperatures, T < T u ,T R ,T d we obtain (see Eq. (|3TJ) ) that 
the free energy has temperature correction proportional to the forth power of the 
temperature ~ T 4 . The high temperature expansion, T ^S> Tu,Tr,T<i given by Eq. 
(|40p reveals the first power of the temperature. Nearby to the sphere the energy has 
the form given by Eq. (|57[) . It has the form of the sum the zero temperature term 
and the Plank form temperature correction. 

The entropy of this system shows interesting behaviour (see Fig. [3]). The entropy 
is positive for any distance between atom and sphere of small radius. It tends to zero 



The thermal Casimir- Polder interaction of an atom with spherical plasma shell 16 

for low temperatures according to the Nernst heat theorem and it tends to constant for 
high temperatures. It has been shown in the Ref. [36 that the entropy is negative for 
short distances between an atom and the flat boundary. To reveal the sign changing we 
found expansion of the entropy over d/R up to the fist power of dj R. The analytical 
and numerical analysis shows that the entropy indeed has the region of the negative 
sign for large radii of the sphere and becomes completely positive beginning with some 
definite value of the d/R (see Fig. |4j . 

As noted above in Introduction the hydrodynamical (plasma) approach has 
limited applications and the more realistic model to describe the interaction between 
an atom with fullerene is the Dirac model. It is based on the 2 + 1 dimensional Dirac 
equation with very small parameter of the mass and the Fermi velocity instead of the 
velocity of light. The future aim of our consideration is the Casimir- Polder interaction 
in the framework of the Dirac model. 
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